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Abstract 

In this paper we clarify the role of heat flux in the hydrodynamic balance 
equations, facilitating the formulation of an Onsager relation within the frame- 
work of this theory. Previously thought to be unobtainable from the present 
form of the theory [X.L. Lei, J. Cai, and L.M. Xie, Phys. Rev. B 38,1529 
(1988)], our verification of the Onsager relation for linear particle and heat 
flux currents driven by electric fields and temperature gradients resolves a 
puzzling issue of long standing. Our results show that, for any temperature, 
when electron density is sufficiently high, the linear predictions of balance 
equation theory exactly satisfy the Onsager relation. The condition of high 
density is consonant with the requirement of strong electron-electron interac- 
tions for the validity of the Lei-Ting balance equations. Our results support 
the validity of this theory for a weakly nonuniform system. We also discuss 
a possible method of extending this theory to a system further removed from 
thermal equilibrium. 
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I. INTRODUCTION 

The balance equation transport theory of Lei and TingS'l was originally developed to 
treat high-field electrical condition in homogeneous semiconductors, and has achieved much 
success in hot-electron semiconductor transport problems. This theory is based on a separa- 
tion of the center of mass of the system from the relative motion of electrons in the presence 
of a uniform electric field. The center of mass is treated as a classical particle, whereas 
the relative electron system, which is composed of a large number of interacting particles, 
is treated fully quantum-mechanically. The theory has been successfully applied to a va- 
riety of transport problems, and the results obtained have exhibited good agreement with 
experiments.! This theory was subsequently generalized to deal with weakly nonuniform, in- 
homogeneous systems by Lei et al.% The resulting hydrodynamic balance equations obtained 
by them consist of the following three equations: (a) continuity equation; (b) momentum 
balance equation; and (c) energy balance equation. 

The form of these hydrodynamic balance equations appears very similar to their classical 
counterparts, generally called hydrodynamic models.i~lii However, in actual fact, they are 
quite different. The latter is derived from the Boltzmann transport equation, as the first 
three moments of that equation. Very recently, the fourth moment was discussed by Anile et 
a/.I^HJ, in an attempt to include the equation describing heat flux. Although, in principle, a 
complete determination of Boltzmann equation is equivalent to the determination of all the 
moments, it is not practical to solve the infinite hierarchy of coupled equations governing the 
various moments. The hydrodynamic approach is based on truncation of this hierarchy af- 
ter the second order moment, and simplification of the remaining equations. However these 
three moment equations by themselves do not form a closed system, requiring input of infor- 
mation about scattering, generally supplied in the form of momentum and energy relaxation 
times. Nevertheless, to accurately evaluate the relaxation times requires a predetermination 
of the distribution function, which is precisely the task that the hydrodynamic models strive 
to avoid. This difficulty is circumvented by one of the following ways. One approach is 
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to calculate the relaxation times by Monte Carlo simulations. Another employs empirical 
forms of relaxation times. The third is to postulate a distribution function with unknown 
parameters, and use the hydrodynamic equations to solve for these parameters. One of the 
most popular parameterized distribution functions is the drifted Maxwell distribution, which 
depends on two unknown parameters, the electron drift velocity and the electron temper- 
ature. The hydrodynamic balance equation approach employs a drifted local equilibrium 
description similar to the latter. In this it employs unknown parameters including the local 
electron temperature T e (R), local electron drift velocity v(R) and local chemical potential 
/i(R). The distinctive features of balance equation theory rest with the ansatz of such local 
equilibrium parameters in an appropriately chosen initial density matrix, which is treated 
quantum mechanically, describing the dynamics of the many-body system of electrons, im- 
purities and phonons. Of course, these unknown parameters are also to be determined from 
the resulting balance equations. It is now believed that the specific quasi-equilibrium form of 
the initial density matrix chosen in balance equation theory is specifically suited to the con- 
dition of strong electron-electron interactions, since it requires rapid thermalization about 
the drifted transport state.00 A salient feature of this new hydrodynamic approach is that 
it includes a microscopic description of scattering in the form of a frictional force function 
due to electron-impurity and electron-phonon scattering, as well as an electron energy loss 
rate function due to electron-phonon interaction. These functions are calculated within the 
model itself, as functions of carrier drift velocity and carrier temperature, along with the 
carrier density, which are themselves determined self-consistently within the same model. 
These hydrodynamic balance equations have recently been applied to device simulations by 
Cai et al.EtB 

A hitherto unresolved question, unanswered since the development of hydrodynamic bal- 
ance equations, concerns the capability of this theory to lead to the correct form of Onsager 
relationsf^H and/or how to determine Onsager relations within the framework of this the- 
ory. There is even some misunderstanding that the energy flux predicted by this theory 
is zero. The purpose of this paper is to clarify the role of heat flux in this theory, and 
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to also show how to generate Onsager relations within the framework of this theory. We 
have closely checked the Onsager relation predicted by this theory and find, that for any 
temperature, when electron density is sufficiently high, the balance equation theory satisfies 
Onsager relations exactly. The condition of high density is consonant with the requirement 
that Lei-Ting balance equations hold only for strong electron-electron interactions. Fur- 
thermore, our results support the validity of this theory in weakly nonuniform systems. To 
our knowledge, this is the first set of hydrodynamic equations which obeys Onsager relation 
exactly. Anile et al. showed very recently^ by Monte Carlo simulation that the Onsager 
relation fails in the traditional hydrodynamic models. 

This paper is organized as follows: In Sec. II we review the derivation of the hydrody- 
namic balance equations. This is not insignificant because we explicitly exhibit the role of 
the energy flux in this theory. Moreover, we also formulate the hydrodynamic force and 
energy balance equations in somewhat different forms than those of Lei et a/.,i which clar- 
ifies the meaning of every term. Then, in Sec. Ill we derive the Onsager relation for linear 
particle and heat flux currents driven by electric field and temperature gradient, and check 
it closely. We present our conclusions and discussions in Sec. IV. 

II. HYDRODYNAMIC BALANCE EQUATIONS 

The starting point of hydrodynamic balance equation theory consists of the following 
fluid-element-composed electron Hamiltonian 



denotes the kinetic energy and Coulomb interaction energy of electrons within a fluid cell 
around R. Macroscopically this cell is small over which all the expectations of physical 
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quantities change little, whereas microscopically it is large enough that a great number of 
particles are within it. and rj are the momentum and coordinate of the i-th electron. 

ifj(R) = XM r <) + - R) (3) 

i 

is the interaction Hamiltonian in which 0(r) denotes the potential of the external electric 
field E, hence E = — V0(r), and $(r) = J2 a u ( r ~ R<0 + J2e u e • V^(r — R^) represents the 
scattering potential due to randomly distributed (R a ) impurities and lattice vibrations (R^ 
stands for the lattice sites). The number density of electrons in the cell around R may be 
written as 

N(R) = £ 5(r, - R) . (4) 

i 

Similarly the R-dependent momentum density is given by 

P(R) = X)Mr i -R). (5) 

i 

Letting v(R) be the average electron velocity in the fluid cell about R, which is a parameter 
to be determined self-consistently from the resulting balance equations, one can write the 
statistical average of the momentum density as 

(P(R)> = mn(R)v(R) , (6) 

with n(r) = (7V(R)), the statistical average of the electron number density. Introducing 
relative electron variables 

p; = Pj -mv(R) , r^ = r,-R, (7) 

which represent the momentum and coordinate of the i-th electron relative to the center of 
mass of the fluid cell around R, we can write the statistical average of H e (R) as 

(# e (R)> = «(R) + ^mn(R)v 2 (R) , (8) 

with 
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«(R) = <£|^«)> (9) 

i 

denoting the average kinetic energy of the relative electron in cell R. It is noted that in 
deriving Eq. (|8]) we have treated electron-electron Coulomb interaction effect in the spirit of 
Landau fermi-liquid theory, which is appropriate for electrons in semiconductors and metals, 
ie., it leads to a self-energy correction in the single electron energy, and also renormalizes the 
bare phonon frequency, jointly with the bare electron-phonon interaction vertex, and also the 
electron-impurity interaction vertexBHil We assume that these renormalized corrections 
are already included in the corresponding quantities. The use of the Hamiltonian above is 
well established and similar to those discussed in the book of Zubarev.il 

Considering the rate of change of particle number density, N(R) = —i[N(TV),H], and 
performing the statistical average, the equation of continuity follows as 

OTl 

_ + V>v) = 0, (10) 

where we have used the relation 

= -i[r h H] =p i / m . (11) 

The particle flux density operator J(R) can be derived from the momentum density operator 
Eq. (D as 

J(R) = ip(R) = £-*(ri-R), (12) 
m m 

and the rate of change of J(R) can be written as 

j(R) = -z[J(R), H] = £ ~(eE + F i )5(r i - R) - V R ■ ]T ^5(r t - R) . (13) 

m ; mm 

Here, we have used the relation 

Pl = H\ = (eE - V$(r<))/m = (eE + F,)/m , (14) 

with Fj representing the force operator of the ith-electron. Transforming to the relative 
coordinate variables, Eq. (0), and performing the statistical average of Eq. (|i~3D , we have 
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|<J(R)> + V ■ «J(R))v) = -V ■ fc**5w) + ^ + ^ • (15) 
OT \ ^ mm/ m m 

where 

<J(R)> = n(R)v(R) , (16) 

and 

f(R) = -(EV^ + R)^)) (17) 

i 

is the frictional force experienced by the fluid cell due to impurity and phonon scattering. 
Since 

(M s ^)=M^ sm ) x= ^ u{r)i ' (18) 

with X as the unit tensor, one follows that 

|<J(R)> + V • «J(R)>v) = -^Vn(R) + ^ . (19) 

This equation can be proved directly to be just the original Euler-type momentum balance 
equation obtained by Lei et alM 

<9v „ 2V« f 

— + v • Vv = — 1 E H , 20 

at 6 mn m mn 

if one takes Eq. ( [TD| ) into account. 

Similarly one can derive the energy balance equation by averaging the Heisenberg equa- 
tion of motion H e (TVj = — i[H e (TV), H], which should be combined with the time derivative 
of Eq. (H) , and yields 

— — h V • \Jh) — _v ' V!i H — mv V • (nv) H — mnv • Vv 
(jij 3 2 2 

-w-v-f. (21) 

Here 

w(R) = W ^ • VStf + R)5(rD) + W V$« + R) • Et*(rJ)) - v(R) • f (R) (22) 

£11/ • lit 
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is the energy transfer rate from electron system to phonon system, and 

J#) = Ef%- R ) ( 23 ) 

^ 2m m 

is the energy flux operator, whose statistical average being 

(Jh(R)> = jj«(R)v(R) + ^mn(R)t; 2 (R)v(R) . (24) 

This is just the energy flux predicted by balance equation theory. Taking this equation into 
account, one can easily recover the original form of the energy balance equation of Ref. f| by 
substituting Eq. (H) into Eq. ([21]): 

— + v-V« = — u(V • v) - w - v • f . (25) 
dt 3 v ; y J 

The resistive force f, the energy transfer rate w, together with the local kinetic energy 
u and the local density n are calculated within the framework of balance equation theory§, 
which requires knowledge of the density matrix p. This density matrix can be determined 
by solving the Liouville equation, idp/dt = [H,p], with an appropriate initial condition. 
In the balance equation theory, the electron-impurity and electron-phonon couplings are 
turned on from t = 0, together with the electric field E. Meanwhile in the present model the 
interactions between different fluid cells are included approximately in the local potential 
with a mean field treatment. Therefore different cells are dynamically independent, and 
thus evolve separately from their own initial state. Thus, the R-dependent initial density 
matrix is chosen such that the relative electron system in the fluid cell is in a local quasi- 
thermal equilibrium state at electron temperature T e (R) and chemical potential /i(R), which 
are parameters to be determined self-consistently from the resulting hydrodynamic balance 
equations; whereas the phonon system is assumed in thermal equilibrium: 

Po = \ exp{-]T[iJ e (R) - v(R) ■ P(R) - p]/T e (R)} exp(-H ph /T) (26) 

Z R 

with H p h and T being the phonon Hamiltonian and temperature. It follows that the resistive 
force and the energy transfer rate are given by 
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f (R) = f (n(R), T e (R), v(R)) = ]T qKq) | 2 n 2 (q, q • v(R)) 

q 

-2^q|M(q,A)| 2 n 2 (q,fi qA -q-v (i ) 

qA 

w(R) =w(n(R),T e (R),v(R)) 

= 2 X; ^qA|M(q, A)| 2 n 2 (q, fi qA - q ■ v d ) 

qA 

with n(x) = (e x — being Bose distribution function; n«, impurity density; f2 q A, the 
phonon frequency of wave q and mode A; tt(q), the electron-impurity interaction potential, 
and M(q, A), the electron-phonon correction matrix element. IT 2 (q, A) denotes the imaginary 
part of electron density- density correction function. Note that f and w depend on R through 
the quantities n(R), T e (R) and u(R). The average local kinetic energy density of the relative 
electrons is 

u = 2 J £e k f[{e^- t j)/T e ] , (29) 

k 

and the local chemical potential /x(R) is related to the local density n(R) of electrons through 
the relation 

n = 2 y £ t f[(e k -ii)/T e ], (30) 

k 

with Ek = k 2 /2m and f(x) = l/(e x + 1) representing the energy dispersion of electrons and 
fermi distribution function respectively. 

There are, altogether, eight variables which need to be determined: the cell velocity 
v(R), the cell electron temperature T e (R), the particle flux (J), the energy flux (Jh), the 
average local kinetic energy density w(R), the local number density of electrons n(R), the 
local chemical potential /x(R), and the total electrical potential 0(R). Moreover, there are 
three balance equations (|T0|), (|19D , and (pl|) , supplemented by four relations (0), (0), (|29|) 
and fl30D, as well as the Poisson equation relating electron density with potential: 

V 2 = -47re[n(R) - n + ] (31) 

with n + as the density of the ionized donor background. These eight equations form a close 
set of equations for the hydrodynamic device modeling. 
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n( 



qA- 



n 



qA 



q-v(R) 



TJR) 



(27) 



,0 



qA 



T 



n 



-^g A - q • v(R) 
TJR) 



(28) 



III. ONSAGER RELATION IN HYDRODYNAMIC BALANCE EQUATION 

APPROACH 

In this section, we demonstrate the Onsager relation, more accurately, we verify the 
validity of hydrodynamic balance equations in regard to the Onsager relation. It is well 
known that the Onsager relation is a manifestation of microscopic irreversibility for any 
statistical system near thermal equilibrium. Therefore any properly formulated statistical 
physics model should satisfy this relation. It is very easy to verify this relation in the 
framework of Kubo linear response theory. Moreover, if one can determine the distribution 
function from the Boltzmann equation, it is also straightforward to verify the Onsager 
relation by calculating the pertinent moments of the distribution function. However, for 
the traditional hydrodynamic verification has been elusive. In fact, in a very 

recent article, Anile et a/.El showed that the Onsager relation breaks down in this model. 
Although they tried to circumvent this difficulty, they did not establish the existence of the 
relation they employed within the model itself by Monte Carlo simulation. Here, we will 
examine the Onsager relation within the framework of the hydrodynamic balance equations. 

The Onsager relationil is concerned with the linear response of the particle current (J) 
and the heat flux (Jq) near thermal equilibrium, which flow as a result of forces Xj on the 
system: 

(J) = L n X! + L 12 X 2 , (32) 
( J Q ) = L 21 X> + L 22 X 2 , (33) 



r 

rl2 r21 



with Xi = — ^V(/i + ecj)) and X 2 = V(l/T). The Onsager relation states that 

(34) 

The heat flux (Jq) relates to the energy flux in Eq. fl24|) through 

(J Q ) = (3 H ) - //(J) . (35) 
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The fluxes (J) and (J#) have already been defined in the previous section by Eqs. (|T6| ) and 



(p4|) . Our first task is to express them in terms of linear response in the form of Eqs. ( |32| ) 
and (H). 

The first relation can be acquired directly by linearization of force balance equation, Eq. 
(|H5|), near thermal equilibrium, so that we only need to consider a steady state with the 
external electric field E and the spatial gradient being very small. Then T e = T and v is 
also very small. We take E, VT and v to be in the x-direction and treat Eq. ( |l9t) to first 
order in the small quantities. This means, for instance, the gradient operator V x = d/dx 
is a first order small quantity and v x is also a first order small quantity, thus V x v x is a 
higher-order small quantity and can be neglected. These facts should be took in mind in all 
of our following calculations. Therefore the force balance equation Eq. (eq2) can be written 
as 



= --?—V x u+ — + — 

Snm m nm 



(36) 



All the quantities in the other two directions are zero. For small v x , f x is proportional to 
v x $ thus proportional to (J x ), and 

fx fx 



p 



n 2 e 2 v x 



ne 2 (J x ) ' 

is the resistivity and independent of v x {(J x )), which is given by 



47T 



n 2 e 2 



E<^|M(q,A)p 



qA 



T 



qA, 



T 



T ' 



g k+q ~ 

T 1 



7tj X ^ I / \ I ^ r \ I 

— L?xKq — n 2 q,w w =o- 

n e q °U) 



(37) 

^(^k+q - + ^qA) 

(38) 



We then have 



/ j \ - — - 2 ^ xU 

\Jx) — o 2 ' 

ep 3 ne z p 

Employing Eqs. ( |29"D and (§t])> we can express Eq. ( |3TS| ) in the form of Eq. fl3"2|), with 



pe z 
T 2 

Ll2 = —2 

pe z 



5^3/2(0 



3*1/2(0 J 



(39) 

(40) 
(41) 
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Here Q — fi/T and the function F u (y) is defined by 

FM = / 7 V : • (42) 
jo exp(x — y) + 1 

The procedure for identifying the linearized heat flux is, of course, similar to that of 
particle flux. Therefore we consider the rate of change of the energy flux operator Jh 
defined by Eq. (|23|) : 

j H (R) = -i[J H (R),H] = -V ■ A 

+ 1 V (eE + FQ- Pi p, _ + 1 ^-(eE + FQp, _ 

2 2m m 2 2m m 

l^pfeE + F,. l^eE + Fip? 

+ - > — -5(ri - R) + - > -— Shi - R) 

2 ^ 2m m 2 4- m 2m v y 

+ 1 V Pi^ + FQ-P,^ _ R) + I V PjP^eE + F,) _ 

2 m 2m 2 m 2m 

where we have used Eqs. flTTp and (|TJ|) again. The tensor A is defined as 

^=£^--%-R)- ( 44 ) 

i 2m m m 



Performing the statistical average of Eq. fl43|) , we have 

5 11 
(3 H ) + V • (A) = (B) + — euE + enE ■ vv + -env 2 E + -v 2 f - wv . (45) 



It is understood that the right hand side of Eq. (|45|) is derived by transforming the coordinate 
and moment operators to the relative variables of Eq. (|7]), before performing the statistical 
averages. The expression of (B) is given in the Appendix, and (A) can be expressed as 

(A) = -(S(R) + uv 2 )l + ( J H )v + v(J# ) - uvv - -mnv 2 vv , (46) 
o z 

with 

s(R) = (E !>.))■ <«) 

This average can be calculated in the balance equation theory mentioned using the density 
matrix p discussed in the previous section, with the result 
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W^E^/(^>. («) 

It should be emphasized here that if the density matrix employed in the balance equa- 
tion is exactly the real physical one, then Eq. ( f45[) should be consistent with Eqs. (fL0D-(plD. 
This is to say that if we have calculated every unknown parameters from the hydrodynamic 
balance equations presented in the previous section, and substitute them in Eq. ([45]), then 
Eq. (fl5|) should merely be an identity. Unfortunately, in actual fact, this is not the case, 
especially when the system is a bit far away from weakly nonuniform system. However, 
here we do not care about it, because we only need this equation holds near thermal equi- 
librium. In this circumstance, the density matrix, chosen in balance equation theory, has 
already been shown to be reasonable, in particular for a system with strong electron-electron 
interactions.0i3 Therefore Eq. (|4"5| ) should yield agreement with the balance equations near 
thermal equilibrium, and we may use it to determine the linear response relation of (Jh) with 
the external forces Xj and examine whether the result obtained satisfies Onsager relation. 

Thus, to the first order in the small quantities, Eq. (|45|) can be written in the form 

^-eu(R)E x - -V^R) + (B x ) = . (49) 
In deriving this equation, we have used the linearized force and energy balance equations, 



Eqs. (^) and fl25|), and (B x ) has also been linearized and is proportional to (J#), which is 
^uv x to first order. Thus we may define 

r n(R)(J H ) ' [M) 

which is also independent of v x ((J#)). Substituting this relation into Eq. (|^) and calcu- 
lating the gradient of S(R) in Eq. (|48|) , we find the average energy flux is given by 



5T 2 F 3/2 (Q T 3 

— -TAj — 

3 m Fi/ 2 {() m 



rX 2 . (51) 



.3F 1/2 (C) 3^1/3(0. 

Subtracting //(J), we obtain the linearized heat flux in terms of Xi and X 2 and can identify 
the linear coefficients of Eq. (|33D as 
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L 21 

L 22 



pe 2 



rpe 2 5F 3/2 (C) 



-T 3 



rn 



m 3Fi /a (C) 
7F 5 / 2 (C) 



-c 



5^3/2(0 



CT 3 



pe z 



5F 3/2 (C) 



3Fi/ 2 (C) 



c 



(52) 
(53) 



.3F 1/2 (C)^ 1/2Uy 3^F 1/2 (C). 

Comparing Eq. (ESI) with Eq. (|4"lD, we find that the condition under which the Onsager 
relation holds is given by 



rpe 
m 



1 



(54) 



We have closely examined Eq. ([54]) for a GaAs system to see if it is indeed satisfied in 
balance equation theory. Both p (Eq. (|38|)) and (B x ) (Appendix) are composed of contri- 
butions due to electron-impurity, electron-LO-phonon, and electron-Ac-phonon scatterings 
(with the electron-acoustic-phonon scatterings due to longitudinal mode acoustic phonons 
via deformation potential and piezoelectric interactions, and transverse mode via piezoelec- 
tric interaction). We have examined each scattering contribution in detail to check Eq. (|54"|) 
separately for each interaction. It is clear that if — e n p ,' / " 1 = 1 is satisfied for each interaction, 



we have 



(l/r)< 

1. Moreover, this procedure effects the fact that the result should be 



independent of impurity concentration and parameters of the electron-phonon interaction 
matrixes. 

The expressions for / obtained from the balance equations are given by 



(55) 



E q s 2 


«(q)l 2 [£ni( q ,a,)]Uo 


(§)(;;) £ q ? 2 


ti(q)| 3 [£n 3 (q, w)] l =0 



due to electron-impurity scattering; and 



'e— ph 



(A) 



E q |M(q, A)| 2 ^ qA (e q + O qA )n^(%^)n 2 (q, fl qA ) 
(|)(^)E q |M(q,A)|^n'(%^)n 2 (q,fi qA ) 
-E q |M(q,A)| 2 ^(^)ni(q,-^ qA ) 
;)(^)E q |M(q,A)| 2 ^'(^)n 2 (q,fi qA ) ' 



+ 



(56) 



due to electron-phonon scattering, for phonons of mode A. I e _ p h(\) is further composed 
of contributions due to electron-LO-phonon scattering, I e -Lo] due to electron-longitudinal 
acoustic phonons by deformation potential coupling, I ed i; and by piezoelectric interaction, 
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I ep i; and due to electron-transverse acoustic phonons by piezoelectric interaction, I ept . TI| 
in Eqs. ( |55| ) and (^) is defined by 



II|(q,a;) = 2vr^£: k (5(£: k+q -£: k + ^) 

k 



. £ k - At \ _ , Z' ^k+q - 



T / \ T 



(57) 



For the LO phonon, Q^lo — = 35.4 meV, and the Frolich matrix element is 
|M(q, LO)\ 2 = e 2 (/t^ 1 — K _1 )f2 /(2£ g 2 ) oc 1/g 2 . (Since the constants in the matrix ele- 
ments cancel in Eq. (|56|), therefore in the following we only specify their relation to q.) 
The matrix element due to longitudinal deformation potential coupling is |M(q, dl)\ 2 oc q, 
that due to longitudinal piezoelectric interaction is |M(q, pl)\ 2 oc (q x q y qz) 2 1 V \ and for 
the two branches of independent transverse piezoelectric interaction: J2j=i,2 \M(<i,ptj)\ 2 (X 
(QxQy + QyQz + QzQx ~ (3q x q y q z ) 2 / q 2 ) / q 5 . For acoustic phonons J7 qA can be written as v s q, with 
the longitudinal sound speed v s being 5.29xl0 3 m/s, and the transverse sound speed being 
2.48 xlO 3 m/s. The effective mass of electron is 0.07m e , with m e denoting the free electron 
mass. 

The results of our numerical calculations are presented in Fig. 1 to Fig. 5, where con- 
tributions to / due to the various interactions discussed above are plotted against electron 
density for several different temperatures. As it is generally believed that the contribution 
of acoustic phonons is important only at low temperature, while the contribution of LO 
phonons is dominant at high temperature, our temperatures are chosen as 10, 20, and 40 K 
for the former, and 50, 300, 500, and 1000 K for the letter. Impurity scattering is present 
at any temperature, so we take T =10, 50, 100, 300, and 1000 K in Fig. 1. From these 
figures it is evident that, for any temperature, when electron density is sufficiently high / 
is exactly unity, indicating that the Onsager relation holds. It is also seen from the figures 
as temperature becomes higher, the electron density needed to make the Onsager relation 
hold is also higher. An interesting exception is the LO phonon in Fig. 2, in which we can 
see that the needed density for T = 300 K is lower than that for T = 50 K, to assure that 

leLO — 1- 
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IV. CONCLUSIONS AND DISCUSSIONS 



In this paper, we have clarified the role of heat flux in hydrodynamic balance equations. 
We have further shown that, for any temperature, when electron density is sufficiently high, 
the hydrodynamic balance equation theory satisfies the Onsager relation. This is consistent 
with the understanding that the Lei-Ting balance equation theory holds only for strong 
electron-electron interactions. Our result supports the validity of this theory in a weakly 
nonuniform system. To our knowledge, this is the first set of hydrodynamic equations which 
satisfies the Onsager relation self contained and without the ad hoc introduction of terms 
which do not originate within the theory. 

However, we should also point out that the hydrodynamic balance equations can only 
be used to describe weakly nonuniform systems. When the temperature gradient is large, 
and/or there is a large heat flux in the system, for example in phenomena as impact ionization 
and heat generation in nonuniform systems, the energy flux equation (Eq. (|45|)), or heat flux 
equation, which is of paramount importance in describing these phenomena, is no longer 
consistent with the other balance equations (Eqs. ([T0|)-([2l|)), and a contradiction emerges. 
This reflects the inadequacy of the assumed initial density matrix, Eq. (p6|), in Lei-Ting 
balance equation theory, by failing to include the detailed information about the physics of 
heat flux. This can be further illustrated as follows: In our deriving the average energy flux 
operator Eq. (p3|) , there should be another term 

on the right hand side of Eq. (|24]). Moreover, in obtaining the average of the tensor A in Eq. 
(H), there should be another term v • (J2i — — m<5( r 0) on the right hand side of Eq. (pEBp . 
These two terms do not vanish when the system is not near thermal equilibrium, and should 
be included in the theory if they are calculated from a real physical density matrix. Anile 
et a/.El have included such terms in their traditional hydrodynamic equations mentioned in 
the introduction. Unfortunately these terms are exactly zero predicted by balance equation 
theory. 
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It is clear that for mediately nonuniform systems and/or systems far from thermal equi- 
librium, an accurate prediction of the behavior of heat flux requires the inclusion of one 
or more additional unknown parameters in the initial density matrix (in high-order terms 
so that they do not violate the particle and momentum balance equations) to be followed 
by their determination from expanded balance equations, which now include the heat flux 
equation(s). This problem is currently under investigation, and the results will be published 
in elsewhere. 
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APPENDIX: 



The expression of (B) is composed of two parts. One is due to collisions with impurities 
((Bj)), and the other is due to interaction with phonons ((B p h)). They are given by 



(Bi) = 2-KUi Kq)| 2 (£ k +q - £k) k4 _ q//2 <Hgk+ q - £ k + q ■ v) 
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(B pfc ) = -4tt £ |M(q, A)| 2 (£ k+q - e k ) k+ J l/2 ^(e k+q - £ k + %x - q • v) 
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FIGURES 

FIG. 1. I due to electron-impurity scattering is plotted as a function of electron density for 
several different temperatures 

FIG. 2. / due to electron-LO-phonon scattering is plotted as a function of electron density for 
several different temperatures 

FIG. 3. / due to electron-longitudinal acoustic-phonon scattering via deformation potential 
coupling is plotted as a function of electron density for several different temperatures 

FIG. 4. / due to electron-longitudinal acoustic-phonon scattering via piezoelectric interaction 
is plotted as a function of electron density for several different temperatures 

FIG. 5. / due to electron-transverse acoustic-phonon scattering via piezoelectric interaction is 
plotted as a function of electron density for several different temperatures 
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